Diffusive energy transport in the 5=1 Haldane chain compound AgVP2S6 
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We present the results of measurements of the thermal conductivity k of the spin 5 = 1 chain 
compound AgVP2S6 in the temperature range between 2 and 300 K and with the heat flow directed 
either along or perpendicular to the chain direction. The analysis of the anisotropy of the heat 
transport allowed for the identiflcation of a small but non-negligible magnon contribution Km along 
the chains, superimposed on the dominant phonon contribution Kph. At temperatures above about 
100 K the energy diffusion constant De{T), calculated from the Km{T) data, exhibits similar features 
as the spin diffusion constant Ds{T), previously measured by NMR. In this regime, the behaviour 
of both transport parameters is consistent with a diffusion process that is caused by interactions 
inherent to one-dimensional 5* = 1 spin systems. 

PACS numbers: 66.70.-(-f, 75.40.Gb 



I. INTRODUCTION 

Anomalous features of transport properties of one- 
dimensional Heisenberg antiferromagnetic (ID HAFM) 
spin S = 1/2 systems were predicted decades ago. In 
particular, it was argued that both the spin and thermal 
conductivity are not expected to be of diffusive char- 
acter as is typical for classical systems without long- 
range order, ^'^ but instead are based on ballistic trans- 
port of spin and energy, at least in the case of ideal 
ID systems. On theoretical grounds it was recently 
demonstrated that this difference is due to the conser- 
vation of spin and energy currents in integrable models 
which apply for ID S=l/2 spin systems. The 
situation is much less clear with respect to the theory 
of transport properties of ID S=l systems. Contrary 
to half-integer spin ID HAFM systems with isotropic 
exchange, for which the spin excitations are S — 1/2 
spinous and the corresponding spectrum is gapless,"'^^ the 
integer-spin chains exhibit gapped excitation spectra and 
5 = 1 magnon excitations. '^'^ Numerical calculations for 
5 = 1, ID HAFM systems relate the gap A to the ex- 
change integral J, such that A = 0.41J.^'* Early theoret- 
ical work predicted that the transport in ideal ID HAFM 
S=l systems must be diffusive.'^ More recent discussions 
of Sachdev and Damle,^^'^^ considering the spin trans- 
port in Heisenberg S=l chains in terms of the nonlinear 
<T- model (NLcrM), also concluded with the claim of diffu- 
sive spin transport. However, Fujimoto,^^ based on the 
integrability of NLcrM, suggested that in a perfect ID 
S=l system, the spin transport is ballistic. The author 
noted, however, that external perturbations, unavoidable 
in real S=l chain compounds and leading to the destruc- 
tion of integrability, may restore diffusive transport. In 
this respect, Fujimoto refers to a previous experimental 
observation of spin diffusion in the S—1 chain compound 
AgVP2S6, by Takigawa et al.,^^ and suggests that the 
spin-phonon interaction might provide this type of per- 



turbation. 

Energy transport via spin excitations, which can be 
probed by measurements of the thermal conductivity, 
has not yet been investigated in S=l chain compounds. 
It seems natural to expect that its features are similar 
to those of spin transport. In a recent paper by Al- 
varez and Gros,^^ relations between the spin and thermal 
conductivities in ID spin systems, similar to the famous 
Wiedemann- Franz law for the ratio between the electrical 
and thermal conductivities in metals, were discussed. In 
the present work, we report the results of an experimen- 
tal investigation of the thermal conductivity of AgVP2S6, 
a compound which is considered as one of the best phys- 
ical realizations of an S* = 1 ID AFM Heisenberg model 
system. AgVP2S6 crystallizes with a monoclinic crystal 
structure of type P2/a.^^ Each ion with spin 5 = 1 
is located inside a distorted octahedron of sulfur ions. 
The zigzag chains formed by these vanadium ions run 
along the a axis. Evidence for an energy gap in the spin 
excitation spectrum was found by measurements of the 
magnetic susceptibility^^'^^'^'^'^'* and inelastic neutron- 
scattering experiments. ^'^'^^ The magnetic properties of 
the compound are well described by the isotropic nearest- 
neighbor Heisenberg Hamiltonian and weak single-ion 
anisotropy 

i? = J^s's^+l + i?^(5D^ (1) 

i i 

with the intrachain exchange constant J/ks = 780 K 
and D/J = 5.8 x 10~^.^^ The interchain interaction 
J' < lO^^J is very weak.^^ The dynamics of the V'^^ 
spins in AgVP2S6 was studied by NMR experiments 
by Takigawa et al.^^ The data are compatible with 
spin diffusion at temperatures exceeding ^ 100 K and 
above a temperature-dependent frequency of the order 
of 10^^ — 10^^ s~^. Using these results, the temperature 
dependence of the spin diffusion constant Ds{T) was es- 
tablished. Below we compare the experimental results 
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for the temperature variation of the spin diffusion con- 
stant Ds(T), probed by NMR, and the energy diffusion 
parameter De{T), which is obtained from the results of 
our thermal conductivity measurements. 

II. SAMPLES AND EXPERIMENT 

For the measurements of the thermal conductivity 
k{T), we grew several single crystals of AgVP2S6 by solid 
state reaction of stoichiometric amounts of Ag, V, P, and 
S, as described in Ref. 20. The largest crystals had di- 
mensions of several mm along the a-axis (the spin chain 
direction), a maximum of 0.5 mm along the 6-axis and, at 
most, 0.2 mm along the c-axis. A single crystalline spec- 
imen with dimensions of 2 x 0.5 x 0.15 mm^ (denoted as 
SI) was used for the thermal conductivity measurements 
along the chain direction. Another sample from the same 
batch served to measure the dc magnetic susceptibility x- 

A standard method for separating the spin-mediated 
thermal conductivity Km{T) from other, less anisotropic, 
contributions to k{T) is based on measurements of this 
quantity both along and perpendicular to the chain di- 
rection. Unfortunately, the dimensions of our crys- 
tals only allowed for measurements of k(T) along the 
chains. For this reason we investigated the heat transport 
also on grain-aligned specimens. Some of these samples 
were previously used for measurements of the magnetic 
susceptibility^^ and for NMR experiments. ^^-^^ Each of 
the specimens with typical sizes of 4 x 2 x 0.2 mm^ were 
single pieces consisting of many equally oriented flake- 
like single crystals. On one such piece (denoted as sam- 
ple Jl), the thermal conductivity was measured along the 
a-axis. Another sample, denoted as J2, with the largest 
spatial extension along the 6-axis, was cut out from the 
central part of sample Jl. Sample J2 was subsequently 
used for experiments with the heat flow perpendicular to 
the chain direction. 

The thermal conductivity was measured using a stan- 
dard imiaxial heat flow method, where the constant heat 
flux along a rectangular-bar shaped sample was produced 
by a heater attached to one end of the sample. The oppo- 
site end of the sample was attached to a copper heat sink. 
The temperature gradient was monitored by a system of 
25 jim diameter Chromcl-Au+0.07%Fc thermocouples. 
We estimate the uncertainty of the absolute value of k 
to be of the order of 10%, because of the uncertainty of 
the sample geometry. The relative error of the measured 
thermal conductivity, however, is only about 0.5%. The 
magnetic susceptibility was measured between 4.5 and 
240 K, employing a commercial SQUID magnetometer. 

III. RESULTS AND DISCUSSION 

The magnetic susceptibility data are shown in Fig. 1. 

The measurements were taken in a magnetic field oi H = 
40 kOe oriented along the a-axis. The data are in good 
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FIG. 1: Magnetic susceptibility x as a function of temper- 
ature. Our results are represented by the solid circles. For 
comparison, previous data reported in Ref. 22 and obtained 
for a powder and a single-crystalline sample are displayed by 
the open diamonds and triangles, respectively. 



agreement with results of Mutka et al.,"^"^ also shown in 
Fig. 1. In particular, the nonlinear X vs. 1/T dependence 
at low temperatures, noticed in Ref. 28 and attributed to 
the weak coupling between flnite-length segments of spin 
chains, is well reproduced in our data (see the inset of 
Fig. 1). The x{T) data can be flt to the equation^^ 

X{T) = ^ exp [-^) + ^ + Xo, (2) 

where the first term on the right-hand side represents the 
contribution from the Haldane-gapped 5=1 chains, the 
second term reflects the above; mentioned flnite segment 
interaction contribution and the last and constant term 
is due to the diamagnctic orbital susceptibility. The best 
fit, shown as the solid line in Fig. 1, is achieved with 
the flt parameters A = 1.77 x 10~^ emu mol~^ K^/^, 
A/ks = 244 K, B = 3.60 x 10"^ emu mor^ K^-^^, 
a — 0.55, and xo = —2.5 x 10~^ emu mol"-'^. The re- 
sulting value for the energy gap A/ks is consistent with 
estimates of 228 K (Ref. 22) and 250 K (Ref. 21) from the 
analyses of magnetic susceptibilities, but somewhat lower 
than the results obtained from NMR^'^ (320 K) and the 
neutron scattering^'"' (300 K) measurements. The power- 
law exponent a = 0.55 is within the range 0.5 < a < 0.8 
established in Ref. 22. 

The temperature dependences of the different thermal 
conductivities are shown in Fig. 2. Each curve exhibits 
a max;imum at temperatures around 10 K and a region 
with a positive slope Ok/OT at temperatures above about 
200 K. The latter feature is presumably due to unac- 
counted heat losses via radiation, typical in standard 
steady-state thermal conductivity experiments probing 
materials with low values of k at high temperatures. Our 
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FIG. 2: Thermal conductivity vs. temperature for two sam- 
ples along the chain direction and one sample in a perpendicu- 
lar direction. The inset demonstrates the slightly irreversible 
behaviour at high temperatures. 



numerical estimates suggest that for our AgVP2S6 sam- 
ples, these effects are negligible only below about 170 to 
200 K and therefore, our data analysis presented below 
is restricted to temperatures below 170 K. 

By rapidly reducing the temperature of the grain- 
aligned samples (Jl and J2) from values in the region 
between 170 and 300 K to temperatures between 80 and 
170 K with a cooling rate of approximately 1 K/s, subse- 
quent slow relaxations of k with time were observed. As 
illustrated in the inset of Fig. 2, even upon a slow vari- 
ation of the temperature (as slow as 0.01 K/s), a small 
hysteresis of k{T) was observed at intermediate tempera- 
tures above the maximum of k(T). We believe that these 
effects are not related to intrinsic properties of the ma- 
terial. Since they were not observed when probing the 
single-crystalline specimen, they are most likely due to 
the multigrain nature of these samples. 

By inspecting the k{T) curves at temperatures below 
60 K, we note that they all exhibit the same qualitative 
features. At higher temperatures, a shoulder- type feature 
in Ka{T), which is absent in the Kb{T) curve, may be iden- 
tified. Similar high-temperature features emerging as a 
shoulder or a second maximum were observed for k(T) 
if measured along the chain or ladder direction of spin- 
chain and spin-ladder compounds. An 
obvious interpretation of this high-temperature feature is 
to ascribe it to a spin- mediated thermal conductivity 
in addition to the common heat conduction via phonons, 
which is reflected in a contribution Kph to the total ther- 
mal conductivity. 

In order to isolate the spin contribution Km[T), accu- 
rate evaluations of other contributions to the total mea- 
sured k{T) need to be made. Since AgVP2S6 is an in- 
sulator, no heat conduction by free charge carriers is ex- 



pected, and only heat transport via lattice excitations 
needs to be considered. Because of the extremely weak 
interaction between spins of neighboring chains, no siz- 
able heat transport via the spin system is expected for di- 
rections perpendicular to the a-axis. Hence only phonons 
contribute to ki,{T), shown in Fig. 2. The shape of the 
Kb{T) curve, typical for phonon heat transport, results 
from the competition between different phonon scatter- 
ing mechanisms, most notably phonon scattering by sam- 
ple or grain boundaries at the lowest temperatures and 
by phonons, defects and magnons at more elevated tem- 
peratures. In contrast to the boundary and defect scat- 
tering which is usually sample-dependent, the scattering 
of phonons by quasiparticlcs in general is expected to 
be sample-independent because it is related to intrinsic 
scattering processes. 

At temperatures much less than A/ks, the number of 
spin excitations that participate in the heat transport re- 
flected in Ka{T) is decreasing exponentially with decreas- 
ing temperature. Therefore, at T <C A/kB the phonon 
heat transport is expected to also dominate Ha{T). While 
the different temperature dependences of k for different 
samples on the left of the k(T) maxima in Fig. 2 can 
be accounted for by differences in the boundary scat- 
tering conditions, it is remarkable that at temperatures 
above the maxima, up to about 40 K, the ratio between 
the K, values for different samples is practically indepen- 
dent of temperature. This implies that the tempera- 
ture dependence of these k{T) curves are identical and 
that in this temperature region, the intrinsic processes of 
phonon scattering dominate. As we have already argued 
in Ref. 29, the anisotropy of the purely phononic heat 
transport is not likely to change at high temperatures. 
This is why we assume that at high temperatures, the 
phonon thermal conductivity along different directions 
and for different samples merely differs by a constant 
factor. With this assumption, we calculated the magnon 
contribution 



K * Kb 



(3) 



with K = 2.7 and K = 2.3 for samples SI and Jl, re- 
spectively. The K values were calculated as the averaged 
ratio Ka/nb in the temperature interval between about 20 
and 35 K. The phonon contributions K * Kb for samples 
SI and Jl are shown by the broken lines in Fig. 2. 

The resulting values of k„j(T) are presented in Fig. 3. 
The uncertainties caused by the subtraction of two ex- 
perimentally measured curves, each with a relative error 
of about 0.5%, are also displayed in Fig. 3. Because of 
the relatively large uncertainties inherent to the evalua- 
tion procedure below about 40 K, it is impossible to draw 
any conclusions about Km{T) in that temperature region. 
Therefore, only data at temperatures above 50 K were 
considered for the following analysis. In the accessible 
regime, Km{T) first increases with increasing temperature 
but tends to pass through a maximum or to reach satu- 
ration at higher temperatures. The maximum absolute 
values of ~ 1 W m~^ for the magnon thermal con- 
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FIG. 3: Magnon thermal conductivity along the a-axis of 
AgVP2S6, extracted from Ka{T) for the two samples SI (a) 
and Jl (b). 



ductivity are one to two orders of magnitude smaller than 
the corresponding values for in S' = 1/2 two- leg spin- 
ladder compounds (La,Sr,Ca)i4Cu2404i'^°'^^''^^ and var- 
ious S = 1/2 spin-chain materials,^®'^^'^^'^^'^'' where the 

enhanced spin-mediated thermal conductivity is thought 
to be the consequence of a quasiballistic energy trans- 
port. 

In view of the implicit interpretation that heat may be 
transported via spin excitations, it seems natural to com- 
pare the energy transport with spin transport. The tem- 
perature dependence of the spin diffusion constant DsiT) 
in AgVP2S6 was measured by Takigawa et al}^ Based on 
our K„i{T) data, we calculated the spin-related energy 
diffusion constant De from De{T) = KmiT) / {C s{T)a'^) , 
where Cs{T) is the specific heat of the spin system and 
a = 2.96 A is the distance between neighboring spins 
along the chains. The low-temperature specific heat of 
Heisenberg type 5 = 1 spin chains was analyzed by Joli- 
coeur and Golinelli^^ employing the quantum nonlinear 
CT-model. The energy s versus the wavevector k, mea- 
sured from the AFM wave vector it /a, is given by the 
equation 



1/2 



(4) 



and the results of numerical calculations at T = are 
Ao = 0.41J and V = 2.49J." The energy diffusion con- 
stant De{T), calculated for the single-crystalline sample 
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FIG. 4: The energy diffusion constant De(T), calculated 
from the thermal conductivity data of sample SI (open cir- 
cles), and the spin diffusion constant Ds{T) from Ref.18 (solid 
circles) are shown in the left panel. In the right panel, the 
corresponding data for the 5 = 1/2 spin chain compound 
Sr2Cu03 are shown. 



SI as outlined above and using Cs(T) as calculated by 
Jolicoeur and Golinelli,^*''^^ is shown in the left panel of 
Fig. 4; the data for the second sample, not shown here, 
arc very similar. The spin diffusion constant Ds{T) from 
Ref. 18 is also displayed in the same figure. The absolute 
values of the two parameters are similar in magnitude at 
high temperatures and they exhibit very similarly shaped 
temperature dependences. This is true although, in com- 
parison with the corresponding feature of DsiT), the 
maximum value of De{T) is about two times larger and 
slightly shifted to lower temperatures. The only other 
spin-chain compound for which data on both the spin 
and the energy transport are available, is the Heisen- 
berg 5 = 1/2 chain compound Sr2Cu03. In the right 
panel of Fig. 4, we present both Ds{T) data from Ref. 40 
and De{T) data calculated from results of our previous 
thermal conductivity measurements on Sr2Cu03.^^ It is 
again obvious that both quantities vary similarly with 
temperature and De ^ 2Ds in the covered temperature 
interval. The only other material for which De and Ds 
have been compared up to now is solid ■^He in the region 
0.05 < T < 0.12 K, where the ratio De/Ds « 2.1 is, 
again, close to two.'*"'^ It is remarkable that indeed the 
factor of 2 difference between De and Ds is expected to 
be a general property of any Bravais lattice with identi- 
cal spins on each lattice point. The fair agreement of 
our results with this expectation thus confirms the re- 
liability of our calculation of the spin-mediated thermal 
conductivity. 

The limited accessible temperature interval for the 
evaluation of Km{T) and the increasing uncertainty with 
decreasing temperatures do not allow a serious analysis of 
its overall temperature dependence. The classical high- 
temperature limit oi De for a spin system is given by^ 
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DEM^J[s{s + i)Y'^/n. 



(5) 



For AgVP2S6, this high- temperature Umit is -D£;,ht ~ 
10"'^^ scc""'^. The observed vahics in the temperature 
range covered in this study are, except at the lowest tem- 
peratures, consistently higher than De,\a, but tend to 
approach this value with increasing temperature. The 
enhancement is, however, much less pronounced than in 
5=1/2 chain and ladder compounds. A better way to 
demonstrate this is by invoking the average mean free 
path of itinerant spin excitations £rn{T), which can be 
defined as £„i = DE/vm, where Vm{T) is the average 
group velocity 



Vm{T) 



de 
dk 



f{e,T)dk 



I 



fie,T)dk 



(6) 



with f{e,T) = (exp(e/fcBT) — 1)"-'^ and the dispersion 
relation e{k) as given by Eq. 4. We calculated i,n{T) 
taking into account the temperature dependence of the 
energy gap, given by the NLcrM at low temperatures as^^ 



A(T) « Ao + (2^AofcBT)i/2 exp(-Ao/fcsr). (7) 

The calculated values of lm{T), shown in the inset of 
Fig. 4, do not exceed 60 A. In contrast, the mean free 
paths of spin excitations in many = 1 /2 chains typi- 
cally reach values of the order of 10^ ^ 29,33,34,35,37 

The data presented above and the subsequent com- 
parison of S* = 1 and 5=1/2 spin chains provoke the 
obvious question concerning the reasons for the much re- 
duced spin-mediated energy transport in the S = 1 sys- 
tem. Two possibilities may be considered. First, extrin- 
sic causes, such as a much stronger scattering of magnons 
by defects and phonons might significantly reduce k„i in 
the compoimd studied here. Second, intrinsic causes of 
reduced energy transport, e.g., as a result of the non- 
integrability of the Hamiltonian even for the pure S = 1 
spin system should be considered. We argue that the 
influence of defects is hardly the major factor leading 
to low values oi De, because the data for the two dif- 
ferent samples of different quality that were studied in 
the course of this work provide very similar results for 
the magnon thermal conductivity (see Fig. 3). That the 
magnon-phonon scattering is the main source for low Hm 
values in AgVP2S6 cannot altogether be ruled out, but 
at present we cannot oSer an obvious reason why this 
scattering should be much stronger in AgVP2S6 than in 
all other spin-chain compounds mentioned above. 

The second and most likely possibility is that the en- 
ergy diffusion observed in our experiments on AgVP2S6 
is intrinsic and is governed by dynamic correlations of 
the S = 1 spin system. Unfortunately, to our knowledge, 
there is no calculation of De{T) available in the litera- 
ture for temperatures much less than J/ks to compare 
our data with. However, a semiclassical consideration of 



5=1 chains by Damle and Sachdcv'^^ and applied to 
AgVP2S6, results in Ds values which are in fair agree- 
ment with the Ds{T) data shown in Fig. 4. Since our 
values for De are compatible with Ds, we take this as 
a clear indication for the same type of behavior of both 
parameters. Nevertheless, the reason for the divergence 
between De{T) and Ds{T) at low temperatures, namely 
below about 100 K is not understood. In Ref. 18, the 
drop of Ds{T) at temperatures below 100 K was consid- 
ered to be the result of a crossover from spin diffusion to 
a region of freely propagating magnons. Our results do 
not contradict this assumption, however a rigorous anal- 
ysis of De(T) in this temperature regime is impossible, 
because of the prohibitively large uncertainty in Km at 
these temperatures. 

A rigorous approach to transport properties of low- 
dimensional quantum spin systems, based on the no- 
tion of a thermal Drude weight Dth{T) in integrable 
spin systems, analogous to the familiar Drude weight in 
the theory of metallic conductivity, has recently been 
developed. i9'^3,44,45,46,47,48,49 rpj^g spin-related thermal 

conductivity Kspin, which here corresponds to a conduc- 
tivity via cither magnons in the 5=1 case (Km) or 
spinous in the 5=1/2 case (k,,), is given as 



*'spi 



n(T) 



lim K, 



spin 



with 



KspinCw) = Dtii{T)6{uj). 



(8) 



(9) 



For an ideal integrable system, Kgpin diverges if Dth is 

nonzero. If perturbations introduce a finite lifetime r, the 
thermal conductivity is also expected to be finite and can 
be represented as Kspin = -DthT- In Ref. 49, where this 
approach was applied to S=l HAFM chains, it is noted 
that this is true only in the case of a large and energy- 
independent r or, equivalently, to a long mean free path 
£. However, impurity scattering, for example, is shown''^ 
to invoke energy-dependent lifetimes, thus leading to a 
more complicated behaviour of k„j which is not simply 
given by a product of the thermal Drude weight and the 
average relaxation time. Since in our case, as reflected 
in the small values of £m which approach the intcrspin 
distance (see the inset of Fig. 4) , the relaxation is quite 
strong, we do not attempt a direct application of the ther- 
mal Drude weight formalism to our result for Km{T) of 
AgVP2S6. This approach seems, however, well justified 
for the analysis of /t^ in 5 = 1/2 ID spin systems, where 
the spinon mean free paths are often rather large. 



IV. SUMMARY AND CONCLUSIONS 

In conclusion, by probing the anisotropy of the thermal 
conductivity of the S=l Haldane-gap HAFM compound 
AgVP2S6, we established the magnitude and the tem- 
perature variation of the spin-related energy transport 
in the temperature region between about 50 and 170 K. 



6 



By comparing the energy and the spin diffusion parame- 
ters De{T) and Ds{T), both derived from experimental 
data, we note that they not only exhibit similar temper- 
ature dependences but that also their absolute values are 
of the same order of magnitude. We argue that the previ- 
ous suggestion of an intrinsic origin of the spin diffusion 
at high temperatures in this system^^ also applies to the 
energy transport that is investigated here. The character 
of the energy transport at low temperatures remains to 



be investigated. 
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